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I. INTRODUCTION 

o ■ 
o : 

While integrating the Einstein-Dirac equations, one has to face two difficulties. 
£SJ ■ The first one is a purely technical difficulty related to the Einstein-Dirac equations being 
,(_» ' a complicated system of nonlinear partial differentiaLeajuations of the second order with 24 
unknown functions. Some particular exact solutionsuEl of the Einstein-Dirac equations in 
homogeneous spaces used to be obtained before only for the diagonal metrics of Riemannian 
, space. 

■ The second difficulty is of a fundamental nature and related to the spinor field in Rie- 
mannian space being definable only in some nonholonomic orthonormal bases (tetrads) to 
be set. In other words, a gauge of the tetrads is necessary. Lots of such gauges are known, 
and different authors proposed different gauges. Taken together these gauges either are 
noninvariant with respect to transformation of variables in the observer's coordinate system 
or are written in the form of differential equations, which results in a complication of the 

, initial system of equations. 

■ Physically, all gauges are equivalent, since the Einstein-Dirac equations are tetrad-choice 
invariant. Mathematically, the use of a bad gauge (i.e. additional equations closing the 

■ Einstein-Dirac ones) may seriously complicate the equations, but the use of a good gauge 
| may significantly simplify them. 

To the-difficulty in a reasonable gauge of tetrads detect is much related that in previous 
paperstrEl solutions of the Einstein-Dirac equations have been obtained only for diagonal 
(3JT), metrics (for such metrics the vectors of the basis of a holonomic coordinate system are 
orthogonal, and hence it is possible to make a natural choice of the tetrads related with an 
orthogonal holonomic basis of Riemannian space) . 

Here we use the tetrad gaugeQ that is algebraic and, at the same time, is formulated in 
an invariant way. Using this gauge allows one to reduce the number of unknown functions 
in the Einstein-Dirac equations by six units conserving the invariance of equations with 
respect to transformation of variables of the observer's coordinate system. 

With the tetrad gauge being used here, all equations are formulated as the first-order 
equations only for two invariants of the spinor field and Ricci rotation coefficients of the 
proper bases defined by the spinor field. The Dirac equations after the tetrad gauge trans- 
form into equations for Ricci rotation coefficients and for invariants of the spinor field, with 
the Ricci rotation coefficients entering in these equations linearly and derivative- free. There- 
fore in a homogeneous Riemannian space the Dirac equations prove to close the Einstein 
equations for the Ricci rotation coefficients without using additional equations. In this case 
one may first integrate the first order equations for the Ricci rotation coefficients and the 
spinor field invariants, and then integrate the first order equations for scale factors. 

With these two circumstances — a reduction of the number of unknown functions by six 
units and a possibility of integrating the second order equations in two steps — considerable 
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simplification of the Einstein-Dirac equations is related, which makes it possible to obtain 
new exact solutions of these equations. 

II. SPINOR FIELDS IN THE FOUR-DIMENSIONAL RIEMANNIAN SPACE 

Let V is the four-dimensional Ricmannian space with the metric signature (+, +, +, — ), 
referred to a system of coordinate with the variables x l and the holonomic vector basis = 
1, 2, 3, 4. The metric tensor of the space V is determined in basis Ji by covariant components 
gij] the connection is determined by the Christoffel symbols T^. Let us introduce in the 
space V a smooth field of the orthonormal bases (tetrads) e a (x z ) (a = 1, 2, 3, 4), by the 
relations 

h aJii Ji — hi G a , (1-1) 

where hi a , h l a are the scale factors. We designate the indices of tensor components in the 
basis Ji by Latin letters i, j, k, ... ; the indices of tensor components in the orthonormal 
basis e a will be designated by the first letters of the Latin alphabet a, b, c, d, e, /. 

The differential of the orthonormal basis vectors e a (x l ) is defined by the Ricci rotation 
coefficients 

de a = A ha b e b dx\ (1.2) 

which are expressed via scale factors 

X.ac = ^[h 3 c {dih ja - djhia) - h 3 a (dih jc - djh ic ) + h^h 1 a h s c (djh sb - d s h jb )]. (1.3) 

Here di = d/dx l is the symbol of a partial derivative with respect to the variable x l . 

Let us determine in the Ricmannian space V a spinor field of the first rank tp{x l ), given 
by contravariant components ip (x l ) (A = 1, 2, 3, 4) in the orthonormal bases e a (x t ). The 
spinor indeces of the components of any spin-tensor are lowered or raised by the formulas 

^ A = e AB i; B , i> A = e AB tp B , (1.4) 

where E = He^sll, E^ 1 = ||e AB || are the covariant and contravariant components of the 
metric spinor, defined by the equations 

7a T - -E la E-\ E T = -E. (1.5) 

Here T is the transposition symbol; -f a are the four-dimensional Dirac matrices that on 
definition satisfy the equation 

7a 76 + 767a = 2 5afci (1-6) 

where I is the unit four-dimensional matrix, \\g a b\\ = diag(l, 1, 1,-1) the covariant compo- 
nents of metric tensor in an orthonormal basis e a . 

It should be noted, that in the equations (1.4) contracting is always done with respect to 
the second index of antisymmetric components of the metric spinor. 

Let us define also the conjugate spinor field with the covariant components = 
by the relation ip + = iji 1 ' f3, where the point above a letter denotes the complex conjugation; 
the invariant spinor of second rank (3 is defined by the equations 

il = -07a/?" 1 , P T = P- (1-7) 

The four-component spinor field ip in general case has two real invariants p, rj, which can 
be defined by the equation 

pexpirj ~ tp + tp + itp + ^ 5 ip, (1.8) 
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where 7 s = jgS abcd j a jbJc7d, e abcd are the components of the four-dimensional Levi-Chivita 
tensor e 1234 = — 1. In the Riemannian space V via the spinor field ip an( A conjugate spinor 
field ip + it is possible to define the proper orthonormal basis e a of the spinor field 

e 1 =n l J t , e 2 =CJi, e 3 = a t J z , e A = u l J t . (1.9) 

The vector components 7r l , £ l , a 1 , u l in the equations (1.9) are defined by the relationsill 



pn l = lm(tp T E^ l ip), 

pC = Re{ip T EY4>), 

pa 1 — ^ + 7*7 5 -0, 
pu 



(1.10) 



where the invariant p is determined by the relation (1.8), the spin-tensors 7* = h l a ^ a satisfy 
the equation 



7V +7-V = 2g lJ I 



(1.11) 



It is obvious, that the scale factors h l a , corresponding to the proper basis e a , are deter- 
mined by matrix 



h\ = 

If the spin-tensors j a , E, (3 are determined by matrices 



tt 1 e a 1 u 1 

n 2 e a 2 u 2 

ir 3 £ 3 a 3 u 3 

4 M 4 4 

7T f <7 7i 



(1.12) 



71 














i 










i 










— i 










-i 






























i 


74 




i 












i 









72 



E = 



1 
0-10 
0-100 
10 

10 
-10 
-1 
10 



73 









i 








- 


— i 











i 








1 











1 


1 








1 









(1.13) 



then the components of |Seinor ip calculated in the proper basis e a , are determined by the 
invariants p, r\ as followal3 



V-' 1= (3, -i/j 2 = i wipexp(|^) 



t/j 3 = 0, 



(1.14) 
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As is well known, the covariant derivative of spinor fields ip, ip + in the Riemannian space 
are defined by the relations^ 



V S V - d s i> - ^As^y V, V s </>+ = d s ip+ + ^+A s ^Yl J - 



For the covariant derivative of spinor fields ip, ip + the equations are valid 



^A S ^ 7 V + hd a Inp- ^ 5 d s rA ip, 
V s ^ + - ^ + ( | A. , + hd s lnp - ^ 5 d sV 



(1.15) 



(1.16) 



Here the invariants of spinor field p, r\ are defined by the equation (1.8), the Ricci rotation 
coefficients Aijk — hj b hk c Ai t b c correspond to the proper bases e a of a spinor field and are 
calculated by the formula 

A s ,ij = 2 (^iV s 7Tj-7TjV s 7ri+£iV s £j-£jV s £i+CTiV^ (1-17) 
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J he. relations (1.16) in the four-dimensional pseudoeuclidean space have been obtained 
The relations (1.16) in a Riemannian space can be obtained from the corresponding 
relations in pseudoeuclidean space by replacement particular derivative to covariant ones. 
The relations (1.16) are satisfied identically in view of definitions of the quantities p, 77, 

A s ^ij • 



III. EINSTEIN— DIRAC EQUATIONS WITH THE COSMOLOGICAL CONSTANT 

Let us consider the set of equations 

7 a V a -0 + imp = 0, 

Rab — ^Rgab + tyab = K-T ab . ^ 

Here ip is a four-component spinor field in the four-dimensional Riemannian space of 
events V, given in the orthonormal basis e a ; m, A, k are constants; R — g ab R a b is the scalar 
curvature of space V, R ab the components of the Ricci tensor calculated in basis e a and 
<7 o6 = diag(l, 1,1,-1). 

The components T ab of the energy-momentum tensor of the spinor field in basis e a are 
defined by the relation 

Tab = 1 (^ + 7aV b ^ - V 6 V + • 7a V> + V> 4 "7b V ' a ip - X7 a i/j+ • 7bV>) ■ (2.2) 

The equations (2.1) are invariant under the arbitrary pseudoorthogonal transformations 
of tetrads e a , therefore for the closure of equations (2.1) it is necessary to add the additional 
equations, defining tetrads e a . Such additional equations usually are called as the gauge 
conditions. Further on we accept the gauge conditions e a = e a , i. e. we accept, that the 
arbitrary tetrad e a in the equations (2.1) coincides with the proper tetrad of the spinor field 
ip. In this case the scale factors h l a in the equations (2.1) coincide with coefficients h l a , 
determined by matrix (1.12). 

The Dirac equations in this case are written as .the equations on the Ricci rotation coef- 
ficients Ka.bc and the invariants of the spinor fielcu 

D a lnp + A b ^ a b = 2ma a sin 77, 
f) a rj + ^e abcd A b , cd = 2ma a cos 77. ^ 

Here D a = h l a di = {■K l di 1 ^ l d il a' l di 1 u l di}- 1 cr a — a a — (0,0,1,0) are the components of 
the vector S3 in the proper basis. 

Coefficients A a , bc m the equations (2.3) are connected with the scale factors h l a by relation 

K,bc = ^ [h j a{D b h JC - bji jb ) + h j c (bah jb + b b h ja ) - h? b (D a h JC + bjlja)} . (2.4) 

The equations (2.3) are the identical record of the Dirac equations in the proper basis e a . 
These equations can be received also from the Dirac equations in (2.1), if to replace in them 
derivatives via the formula (1-16) and to produce algebraic transformations. 

The replacement in equations ^2.3) coefficients A a . bc via h l a yields the following system 
of the invariant tensor equationsQ 

Vj/97r* = Q, Vip^ = 0, V l pa i = 2mpsinr], V i pu i = Q, 

V *?7 - ^S V " lS {■K.jVmlTs + £jV m £,s + (TjV m (T s - UjV ' m u s ) = 2ma l COS 77. 

The Einstein equations in the proper basis e a it is convenient to write as 

Rab = nf ab + (^nmpcosri + A j g ab . (2.6) 
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For transformation of the Einstein equations to the form (2.6) it is necessary to take into 
account, that in view of the equations (2.1) the follow equation is valid 

R = -nT a a +4A = Kmpcos?7 + 4A. (2.7) 

The expression of tetrad components of the energy-momentum tensor in gauge e a = e a 



is 



obtained i 



rp 1 



-v b D a r) - a a b b rj + i<r e (^A a . cd e b cde + A b , cd e a ' 



de 



The tensor components R ab we can express via the Ricci rotation coefficients 
1 



R ab 



--0, 



(/>'•, A,.,/ 



hhA r 



V,6 



'A r 



f 



' : A 



f,b 



(2. 



(2.9) 



Here g = det \\g i3 \\. 

Let us give also the expression of the Ricci rotation coefficients A Qj t> c immediately via the 
vector components iti, £i, cr.;, Ui, which is obtained from the formula (2.4), if replace the 
components h % a by the formula (1-12) 



A M2 

A3,12 
A4 ; 12 
Al,23 
A2,23 
A4,23 
Al,31 
A2,31 
A 4 ,31 
Ai 4 4 
A2,14 
A3 J 14 
Ai ; 24 
A2,24 
A3,24 
Ai, 3 4 
A2,34 
A3 : 34 



CDnn - nW 2 Tn, A 2)12 = -tt 1 D 2 ^ + e^e,, 

|(e£» 4 ^ - nWit-) + K-vtb^i - nw 2 u t + + eb lUi ), 

§(o\Di6 - CDxai) + K^Da^ + aW^ - tvW^ - C^tt*), 
^D^i ~ ebati, A 3 ,2 3 = -CD 3 ai + a l D 2 a u 
\{a l D^ - eb 4 <r t ) + \{a l b 2Ul + b 2 (Ti - ?b 3Ui - uW 3 &), 
-a l Di-Ki + TT l b 3 TTi, A 3j 3i = ir l b 3 ai - o % b 1 o l , 

±(^£ 2 ^ - <J l b 2 ^) + l^b^i + cb^i - zbxai - aw^), 

\{jr l b 4 ai — o % b±-Ki) + ^(j: l b 3 Ui + u l b 3 -Ki — a % b\Ui — u l b\Oi), 

U l D\-Ki — ^D^TTi, A4 ; l4 = —^D^Ui + U l DlUi, 

K^DaTTi - ir l D 2 u t ) + \{CD lUl + u l b^ - ir'D^ - fL^), 
\{u % b 3 iti - ir l b 3 Ui) + \ (a l DiUi + uWia, - -K l D 4 ai - aW^i), 
^(rjb^i - CDxUi) + \{^b 2Ul + uW^i - ttW^ - fD^), 

u l b 2 ^ - A 4 ,24 = -eb 4Ui + uw 2Ul , 

|(v*& 8 6 - Cb 3 u t ) + ±(aW 2 u t + uW^ - CD 4 ai - <rW 4 &), 
Tifabxdi — a l biu t ) + ^(TT l b 3 Ui + u l b 3 -K l — irW^ai — a l b i -K i ), 

\{u i b 2 u l - a l b 2Ul ) + ±(eb 3 ui + U w 3 ^ - eb 4 ai - <r i b 4 &), 

u l D 3 ai — a l D 4 ai, A4 34 = —a l D 4 Ui + u l D 3 Ui. 



(2.10) 



At an integration of the Einstein-Dirac equations it is useful to take into account also the 
matrix of the components of the energy-momentum tensor, which is obtained in correspon- 
dence with definition (2.8): 



2Ai j2 4 A 2j2 4 — Ai,i4 —DlT] + A 3j2 4 A 4j2 4 — Ai j i2 

A2,24 — Ai j i4 — 2A 2: 14 —D 2 rj — A344 — A 2 ,12 — A4 ; i4 

-Z?l77 + A3 i2 4 -■p2?7-A3 ! i4 ^ 2 ^3^ -Z?47?-A 3; i2 

A 4: 24 — Ai 12 — A242 — A4^i4 —D 4 T) — A342 — 2A4 : i2 



(2.11) 



The equations (2.3), (2.6), (2.9) — (2.11) in the given system of coordinate x 1 represent 
a closed set of equations for determining the functions TTi(x^), £i(x J ), Ui{x^), Ui(x 3 ), p(x 3 ), 
r](x J ). The contravariant components of the metric tensor of the Riemannian space are 
expressed via ir^x 3 ), ^{x^), (Xi(x J ), Ui(x J ) by the relation 



9 lJ - h\hhg ab 



a a J 



u u J 



(2.12) 
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IV. GENERAL EXACT SOLUTION OF THE EINSTEIN— DIRAC EQUATIONS 

IN HOMOGENEOUS SPACE 



Let us consider a four-dimensional Ricmannian space referred to synchronous system of 
coordinates with variables x l , in which on definition are valid the equations 

544 = .9 44 = -l, 54a =. 9 4 « = 0, a = 1,2, 3. (3.1) 

We shall seek the solutions of the equations (2.3), (2.5), (2.6), (2.9) — (2.11) in the syn- 
chronous system of coordinate in the supposition, that all unknown functions depend only 
from parameter x 4 = t. Thus, it is suposed, that the space V is the homogeneous Bianchi 
1 type Riemannian space. 

In this case the equations (2.5) are written as 



^(V^Ptt 4 ) = 0, di{^ P e) = 0, d4V^9P^) = 0, 
di{\f~9 P® 4 ) — 2mv / — 5/9 sin 7/, d^rj = — 2ma 4 cosr/. 



(3.2) 



The set of equations (3.2) and equation 

44 4 4 t4|4 i JJ ,4.4 (o o\ 

are the complete set for determining the quantities it 4 , £ 4 , a 4 , u 4 , r\, p\J—g. The general 
solution of the equations (3.2), (3.3) has the form 

C p 7T 4 _ U 4 1 



PV^9 C v C s C u ^l + C^cos 2 (2mt + ^y , g 4) 

eC a sin(2mt + <p) . 1 + iC a cos(2mt + <p) \ • ) 



y/l + C 2 cos 2 (2mt + (pY £ ~^l + Clcos 2 {2mt + <py 

where ip, C n , C|, C a , C u > 1, C p > are the integration constants; coefficient e can take 
any of the two values +1 or —1. By virtue of the condition of synchronism (3.3) the constans 
C are connected by the relation 

(CV) 2 + (Q) 2 + (C ff ) 2 - (C„) 2 = -1. (3.5) 
Let us introduce the notation h a — h 4 a — (it 4 , £ 4 , ct 4 , u 4 ). From the solution (3.4) we find 

{hlMM}= 1 ={C«,Ct,C u }. (3.6) 

y/l + C£ cos z (2mt + ip) 

It is clear from here, that by virtue of solution (3.4) the directions of the three-dimensional 
vector with components hi, h%, hi does not depend on parameter t. 

From definition (2.4) it follows, that the Ricci rotation coefficients A a M can be presented 
in the form 



where by definition 



A a , 6c = ^(h b s ac - h c s ab - h a a bc ), (3.7) 



s ab — s ba — h % ad^hn, + h l b dihi a , , 

U . V W . W [ O.C 

a a b — -aba = h l a dihib - h l bdih la . 



The quantities h a in the formula (3.7) are determined by solution (3.4) as the functions of 
parameter t, therefore the relation (3.8) expresses 24 dependent functions K a ^ c only through 
16 functions s ab , a ab . 

By virtue of definition (3.8) the quantities s a b in a synchronous system of the coordinates 
satisfy the identities 

h b s ab = 0, s a a = 28 A \n^g-. (3.9) 
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From the first identity in (3.9) it follows, that the ten symmetrical components s ab contains 
generally the six independent ones. 

From the equations (2.3) it follows, that the antisymmetric quantities a ab are determined 
by equality 



i ab = 4m [(a a h b - a b h a ) sin rj - s abcd ^ c h d cos 77] 



(3.10) 



Using the equations (3.9), (3.10) and definitions (2.9), (2.11) for T ab , R ab we can transform 
the Einstein equations (2.6) into the equivalent set of equations 

dii^f^gsah) - 2^g (tocos ry+ \np) (e ce f a s b f + e ce f b sj) h c a e - 
-2niy^g(h a s bc + h b s ac )a c sinr) = (urn py/^g cost) + 2\^/~ r g){g ab + h a h b ) , (3.11) 
(s a a ) 2 — s ab s ab — 8(Kmpcosi] + A). 

The quantity p^/—gcosrj in the right part of the first equation (3.11) by virtue of the 
solution (3.4) is constant 



P\ 



-g cos rj — eC p 



(3.12) 



The first equation in (3.11) is obtained by contracting of the Einstein equations in (2.6) with 
the tensor components 6® + h c h a by the index a. The second equation in (3.11) is obtained 
by contracting of the Einstein equations in (2.6) with the components of tensor g ab + 2h a h b 
by the indices a, b. 

The contracting of the equations (3.11) with g ab by the indices a, b gives the equation 



8484 



-<J 



-nmeCp + 3A\/— g 



(3.13) 



for definition of the quantity y/—g. 

If the cosmological constant A is positive A > 0, the solution of the equation (3.13) has 
the form 



— sum ^ 



1 



fi sinh(V3At) + / 2 cosh(v / 3At) 



(3.14) 



where fx, /a are arbitrary constants. 
If A < 0, for y/—g is obtained 



-9 = 



envn 
"2A~ 



Cp{-l + /ain[V=3A(t-t )]} 



(3.15) 



where /, to are the integration constants 

The case A = has been considered inlii 

As a direction of the three-dimensional vector with components h\, I12, /14 does not depend 
on parameter t, then by the constant Lorentz transformation of the basis vectors ei, 62, 64 
it is possible to transform the components hi, /i2 to zero. The initial set of equations (2.3), 
(2.6), (2.9) - (2.11) is invariant under the arbitrary constant Lorentz transformation of the 
vectors ei, 62, 64. Therefore without loss of generallity it is enough to consider a solution 
of the equations (2.3), (2.6), (2.9) - (2.11) only at hi = h% = 0. Under this condition the 
first equation in (3.11) can be written in the form 
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d4(s 33 ^g) - 4ma 4 sin 77 s 33 ^/=g = (enmC p + 2A< v /=#)(? 



4\2 



di(s2 3 \/—g) — 2ma 4 sin 77 s 23 \/~ r g — ( 2m cos 77 + \np ) w 4 si3\/— g = 0, 

c>4(si3\/— ) — 2mcr 4 sin 77 s 13^/^9 + I 2m cos 77 + i/cp 1 7i 4 S23-\/-j? = 0, 

cMsn\/— 5) + 21 2m cos 77 + j/tp )w 4 si2\/— 5 = enmC p + 2A v ^ z g, 

di{s22\f—g) — 21 2m cos 77 + i«/7 Ju 4 Si2\/— 5 = enmC p + 2A v ^ z g, 

d i {si2^/ ZI 9) - ^2mcos77 + j/tpju^Sn - s 22 )\/ = 5 = 0, 

di{sny/—g ) + I 2m cos 77 + ^/tp )u 4 S24\/~ 5 — 2mu 4 sin 77 S31 5 = 0, 

^4(824^/— ff) — ( 2m cos 77 + ^/tp ]7i 4 si4\/^5 — 2mu 4 sin 77 s 2 3\/— 5 = 0, 
dn(s 3 4^J ~ g ) — 2mcr 4 8^77534^—3 = u 4 (sKmC p a 4 + 2m sin 77333^/— g), 
<94(s44-\/ :r g) — 4mu 4 sin 77 534-^/^ = (enmC p + 2X^/^g) (cr 4 ) 2 . 



(3.16) 



The general solution for s a b has the form: 



sn = pu* 

S22 = pu 4 
2 



3 3CpC u 

1 „ 2 



-.g+-Bsin2(C + /3) 



3 3 Cp Cu 



p(u 



4\3 



N 



S33 

s i4 = ^pu\a 4 f[N + 



C P C U 



d 4 V^--Bsin2(C + /?) 



C P C U 



S12 = --pu 4 Bcos2{( + f3), 



S3i =~-pcj\u 4 f\N + 



C ' P C U 



(3.17) 



3 f 

S13 = -p(u 4 ) 2 Acos(C + a), S14 = -pu 4 cr 4 Acos(C + a), 

I 1 

S23 = — p(w 4 ) 2 Asin(C + a), s 24 = -p7i 4 cr 4 yl sin(C + a), 

where A, B, TV, a, (3 are the arbitrary constants; the quantity \J—g in the equations (3.17) 
is defined by the equalities (3.14), (3.15); the quantity ( is defined by the follow relation 



2mcos?7 + ^-np ) u'dl = rarctg 



(tan(2m7; + ip) 



1 



(3.18) 



The integral r = J pu 4 dt in the formula (3.18) depends from the value of cosmologic constant 
A and will be calculated further. 

Substituting in the second equation in (3.11) the components s a b by the formulas (3.17), 
we obtain a connection between the integration constants A, B, N and /1, f2- In the case 
A > is obtained the following relation 



fi fi = 1 



A 



1 



1 



1 



Ct -A + -B z + -N z < 1. 



It is obvious, that if /| > / 2 , then the formula (3.14) for v / ~5 can be presented as 

, £K77l 



2A 



; C p {-l + /cosh[V3A(t-7; )]}, |/|<1, 



(3.19a) 
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where t , f are the arbitrary constants, f 2 = /f — j\ . 

The singular points of the solution (3.19 a) are determined by the equation / cosh [\/3A (t— 
t )] = 1. The solution has one singular point, if / = 1 and two singular points, if < / < 1. 
If / < (in this case e = —1), the solution has no singular points. 

If fi < fit f° r V~~9 the equality is valid 



'-9 = 



enm 



2A 



C p {-l + /sinh[\/3A(i-t )]}, 



(3.196) 



where f 2 = ff — f 2 . In this case the solution has one singular point, defined by the equation 
fsinh[V3X(t-t )] = 1. 

If A < 0, the following formula for the integration constant / in (3.15) is obtained 



d LA 2 + -B 2 + ^N 2 )> 1. 



/ 2 = 1 



In this case there are indefinitely many singular points, are determined by the equation 

The condition \J~ r g > yields a restriction on the values of the integration constants and 
the domain of existence of the solution. 

Using definition (1.14) and formulas (3.4), we can write out solution for the components 
of spinor field in the proper basis 





i> = ± 



1 + iC a cos(2mt + tp) 



he. 



o 



1 — iCcr cos(2mt + <p) 



IsC. 
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where y/—g, depending on a sign A, is defined by the formula (3.15) or formulas (3.19 a,b). 

The equations (3.4), (3.10), (3.17) completely determine the Ricci rotation coefficients 
A a ,bc by the formula (3.7) and represent the first integral of the equations (2.3), (2.5), (2.6), 
(2.9) — (2.11). For determining of the general solution of the equations (2.3), (2.5), (2.6), 
(2.9) — (2.11) it is now enough to integrate the equations (3.8), from which follows 



1 

d 4 hi b = -(s ab + a ab )hi a . 



(3.20) 



In view of definition (1.12) the quantities hi a and found solutions (3.10), (3.17) for s a b, a abl 
the equation (3.20) is possible to transform to the form 

■^j- (u 4 (Tj — o^Uj) — -^A\^j cos(C + a) + £j sin(C + a)] , 

Id, 1 



+ (u 4 aj — a 4 



3y/—g dr 

d , 1 2m 

- i-V~9 ~7i * cos T] 

g dr 6 p 



(3.21) 



- ^fe - i7T J -) J Bexp[-2i(C + (3)] + l -A{u 4 a - a 4 u 3 ) cxp[-i(C + a)]. 

The differential dr in (3.21) is defined by equality dr = pu 4 dt. The set of equations (3.21) 
should be supplemented by the condition of synchronism 



g>4 a = <Ji<j a — UiU a = 0, a = 1, 2, 3. 
For the parameter r we have 



(3.22) 



J pu 4 dt = C P C U J ■ 



dt 
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If the cosmological constant is positive A > and Km ^ 0, a calculation by means of the 
equations (3.19 a,b) gives 



2XC U 



2 Vac 



/- 



dt 



1 



and 



T = 



sum yj\ + f 2 
2XC U 



fsmh[V3X(t-t ) 
In 



/ - VP + 1 + tanh \ V3A (t - to) 



f + x/FTT + tanh (t - t ) 



(3.23a) 



to 



dt 



sum J -l + /cosh[\/3A(i-i ) 



2V\C U 



1 



In 



/ - 1 + \/l - f 2 tanh 5^3A (t - t ) 



/ — 1 — Vl-/ 2 tanh±\/3A(t-t ) 



(3.236) 



+ T , 



where To is the arbitrary constant. 

At A < (f 2 > 1) by means of (3.15) we obtain 



T = 



2XC U 



sum 



I- 



dt 



+ /sin[V=3A(t-io) 



2 V A C u 



enm 



In 



tan i^3A (t -t )-f- x/f 2 ^! 



tan i^3A (i - i ) - / + v 7 /^ 



(3.24) 



+ T . 



Let us change now the unknowns functions (tt\, £\, a x , u x ) — > (7r", <?a) m the equations 
(3.21) 

^A+iTTA = (^+^A)(V^) 1/3 CX P (-iAT-zc), 

^ = e A (V^) 1/ Vexp('-iivA (3.25) 

Here £ is defined by the relation (3.18), quantities er 4 , u 4 are defined by the solution (3.4), 
\J—g is determined by equalities (3.15) or (3.19 a, b). In outcome the condition of synchro- 
nism (3.22) is satisfied identically, and the set of equations (3.21) goes over into the set of 
linear equations with the constant coefficients: 



-2i/3 



^-6 X = ^A(TT° x cosa + f x sma)+^N6 x . 



-A6 x e 



(3.26) 



At j = 4 the equations (3.21) are satisfied identically by virtue of conditions hi — h 2 = 0. 
Using definition (2.12), we obtain the expression for the spatial components of metric 
tensor g a p via the vector components 7r°, 



9*0 = ( xf~9 ) 3 e-i ^ (tt° t$ + C § + a e p ) . 



(3.27) 



To the equations (3.26) for each value of an index A = 1, 2, 3 corresponds the characteristic 
equation for an eigenvalue q 



det 



B sin 2/3 -B cos 2/3 + k A cos a 
-B cos 2/3 — k —B sin 2/3 A sin a 
A cos a A sin a 2N 



ql}=0, 



which is possible to write as 

2(AT - 2q)(16q 2 + k 2 - A 2 - B 2 ) + A 2 [2N - Bsin2(a - j3)} = 0. 



(3.28) 
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The solution of the equation (3.28) in general case is determined by the Cardano formulas. 
A simple solution of this equation is obtained in the case A = and 2Ai,p^ B sin 2 (a — 0). 
A solution of equations (3.26) in the case A = has been considered int-lo. 

Let us consider the case, when is valid the equation 2N = B sin 2(a-/3),i^0. In this 

1/2 

case for the eigenvalues we have q\ = ^N, q 2 — \{A 2 + B 2 — ft 2 ) and 93 = — j(A 2 + 
1/2 

B 2 — K 2 ) . If < A 2 + B 2 < ft 2 , it is possible to write a solution of the equations (3.26) 
as follows 



t\ + i7T l = e- ta {-AQ x ei Nr + [-{k + Be 2l{a -^)F X +4iAG A ]cosAr 

+ [—(ft + Be 2i ( a ~^)G\ — AiAF\] sin At}, (3.29) 
X = [k — £> cos2(a — f3)\Q\e^ NT + A{F\ cos Ar + Gx sin Ar). 



Here A — k 2 — A 2 — B 2 ; Fx, Gx, Qx are the integration constants connected by the 
relation 

e al3x F a G p Qx = (k 2 -A 2 - B 2 )~ 3 / 2 , (3.30) 

where e a ^ x is the Levi-Chivita symbol. A relation (3.30) expresses the equality between the 
value of y/^g, calculated from the solution (3.29), by its values (3.15), (3.19 a,b). For spatial 
components of the metric tensor by means of the formulas (3.27) the following expression is 
obtained 

g a p = (V 3 5) 2/3 e-^ T {Q Q Q^{v4 2 + [ft-Bcos2(a-/3)] 2 } e ^ 

+ (F a Gp + FpG a ) [M sin 2Ar — 4AB sin 2(a — 0) cos 2Ar] 

+ ^ Q J>[16A 2 + Af + Mcos2AT + 4A J Bsin2(a-/3)sin2Ar] 

+ G a G p[16A 2 + M - Af cos2AT-4ASsin2(a-/3)sin2AT] { - 3 ' 31 ' 

+ (F a Q p + F p Q a )2KAei Nt cos At 

+ (G a Qp + G p Q a )2 K Aei NT smATy 

where the contracted notation is entered 

M = A 2 + B 2 + kB cos 2(a - /3). 

The quantity y/—g in the right part of equality (3.31) is defined by the equations (3.19 
a,b) or (3.15). It is obvious, that by the constant transformation of the variable ar, x 2 , x 3 
it is always possible to transform components F a , G a , Q a to the form 



F a = {{k 2 -A 2 -B 2 )- 1 / 2 ,0,0}, G a = {0,(k 2 -A 2 -B 2 )- 1 / 2 ,0}, 
Q a = {0,0,(K 2 -A 2 ~B 2 )-i/ 2 }. 



(3.32) 



The formula (3.31) at A 2 + B 2 ^=0 determine the oscillatory regime of an motion to singular 
points of the solution. 

If A 2 +B 2 > k 2 , then the trigonometric functions in (3.29) are replaced by the hyperbolic 
ones. 

To the diagonal metric 5^ there corresponds the case, when the integration constants A, 
B in the equations (3.26) are equal to zero. In this case we have 

£a + i7r l = K h F \ + iGx) cxp \kt, , 3 



= kQx exp o Nt 



Taking into account equality (3.27), we obtain the solution for the components of the 
metric tensor 



2, , \2/3 

g a p = « (y-g) 



" NT (F a Fp + G a Gf3)+e3 NT Q a Qp . (3.34) 
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Here the quantity \J—g is defined by the formulas (3.15) or (3.19 a,b), and e Nr is defined 
by the relations 



2f 



eNT= f-i + VT^p^HV^it-to), for A>0; n<fl 



f-1- V 1 - f 2 tanh ±V3A(i-t ) 



2 s 



eNT= f-VPTl + ^V3X { t-t ) for A> 

VZ + v/T^+l + tanhiVSA^-to)/ 



, r _ . tan§v/=3A(t-t )-/- y/f*=l 



2 s 



e NT = -»> _L v ; for A < 0. 

^taa|V=3A(* - *o) - f + ^P^T J 

If Fa, G\, Q\ to define by the relations 

F a = {kT 1 , 0,0,}, G a = {0,K-\0}, Qa^M,* -1 }, (3.36) 

then the metric (3.33) is diagonal. 
J£A = B = N = 0, then 

gap = K 2 {V=g) 2/3 {F a F p + GaGp + Q a Qp) , (3.37) 

where \J—g is determined by the formulas (3.19 a), (3.15) and f 2 = 1. 

Some particular solutions of the Einstein-Dirac equations with the cosmological constant 
in the homogeneaus Riemannian space of the Bianchi Ltjfpe-sjith the diagonal metric have 
been obtained inErQ (without cosmological constant seeuQ'olij) . 

Let us consider the transformation of coordinates (x a , t) — > (x a , r) defined by the equation 
dr = pu 4 dt. The function r(i) is determined by relations (3.23 a,b), (3.24). Quantities 
s/—g g Al in the transformed system of coordinate with variables x a , r have the form 

,4a\/ _ n I i — ~ „44\/ 



{V~99 Aa )' = 0, (V^55 44 )' = const. 

From here follows, that in the system of coordinate with the variables x a , r is satisfied the 
harmonicities condition 

ft(V^')] / = -f(V=s3 4i ) , = o. 



dr 

Thus, parameter r is the time in the harmonic system of coordinate. 
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